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基于 GASA-FCM 混合 聚 类 与 霍 夫 变换 的 欠 定 混合 矩阵 估计 ， 


李 虎 , 徐 岩 
(兰州 交通 大 学 电子 与 信息 工程 学 院 ,， 兰州 730070) 


摘 要 : 针对 模糊 C- 均 值 聚 类 (fazzy C-means clustering，FCM) 算 法 在 欠 定 混合 算 阵 估计 中 精度 低 、 鲁 棒 性 差 的 缺点 ， 
提出 一 种 基于 遗传 模拟 退火 优化 FCM(GASA-FCM) 混 合 聚 类 和 霍 夫 变 换 的 欠 定 混合 矩阵 估计 算法 。 该 算法 首先 结合 了 
模拟 退火 算法 (Simulated annealing algorithm，SA) 全 局 搜索 、 高 精度 的 优点 和 遗传 算法 (genetic algorithm，GA) 强 大 的 空 
间 搜 索 能 力 ， 将 经 遗传 模拟 退火 算法 得 到 的 聚 类 中 心 点 赋 给 FCM， 避 和 免 了 初 值 选 择 的 随机 性 。 再 利用 霍 夫 变换 对 聚 类 
得 到 的 每 一 类 数据 的 中 心 进行 修正 ， 提 高 混合 和 琵 阵 的 估计 精度 。 实 验 结果 表明 ， 提 出 的 算法 明显 改善 了 算法 的 稳定 性 
和 混合 矩阵 估计 精度 ， 具 有 一 定 的 有 效 性 和 可 行 性 。 
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中 图 分 类 号 : TN911.7 doi: 10.3969/j.issn.1001-3695.2017.10.0932 


> Underdetermined mixing matrix estimation algorithm using GASA-FCM based mixing 
(CO clustering and Hough transform 


Li Hu, Xu Yan 
(School of Electronic & Information Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China) 


Abstract: With regards to the low precision and poor robustness of fuzzy C-means clustering algorithm (FCM) in the 


underdetermined mixing matrix estimation, this paper proposed an algorithm using genetic simulated annealing optimized 


FCM(GASA-FCM) based mixed clustering ”with hough transform to solve underdetermined mixing matrix estimation. It 


could combine the global search, high-precision advantages of simulated annealing algorithm(SA) and powerful search ability 
of space of genetic algorithm(GA), FCM assigned the clustering center point obtained by genetic Simulated annealing algorithm, 
which avoided the randomness of initial selection. The center of each kind of data obtained by clustering by hough transform is 
modified to improve the estimation accuracy of the mixing matrix. Experimental results show that the proposed algorithm 
significantly improves the stability of the algorithm and the accuracy of the mixing matrix estimation, and has certain validity 
and feasibility. 

Key Words: fuzzy C- means clustering algorithm; genetic algorithm; simulated annealing algorithm; Hough transform; mixing 


matrix estimation 


矩阵 估计 的 稳定 性 和 精度 都 不 高 户 。 针 对 FCM 聚 类 算法 存在 的 


0 引言 
问题 ， 文 献 [3] 提出 基于 遗传 算法 优化 的 FCM， 利 用 遗传 算法 
基于 信号 稀 玻 分 析 理 论 提 出 “两 步 法 ”的 欠 定 盲 源 分 离 的 全 局 搜索 性 能 ,， 求 得 初始 聚 类 中 心 , 避免 了 FCM 由 于 聚 类 中 


法 中 ， 解 决 问题 的 核心 是 混合 矩阵 的 估计 ， 因 为 混合 矩阵 估计 心 选择 不 当 对 聚 类 结果 的 影响 。 文献 [4] 提出 差分 进化 的 方法 来 


的 准确 性 决定 源 信号 分 离 的 精确 度 。 对 于 欠 定 混合 矩阵 估计 ， 优化 模糊 C- 均 值 聚 类 ,避免 了 FCM 依赖 初始 值 ， 容 易 陷入 局 部 
研究 者 利用 观测 数据 在 时 域 或 时 频 域 呈现 的 直线 聚 类 特性 ， 提 最 优 的 缺陷 。 

出 了 一 系列 的 聚 类 算法 ， 通 过 聚 类 结果 计算 直线 向 量 ， 对 混合 本 文 提出 基于 GASA-FCM 混合 聚 类 与 霍 夫 变换 的 欠 定 混 
和 矩阵 进行 估计 吕 。 目 前 运用 最 广泛 的 是 模糊 C 均值 聚 类 算法 ， 合算 阵 估计 算法 ， 在 把 握 全 局 搜索 能 力 比较 强 的 遗传 算法 中 融 


它 是 硬性 C- 均 值 的 自然 推广 ， 计 算 简单 方便 ， 易 于 实现 。 由 于 入 局 部 搜索 能 力 较 强 的 模拟 退火 算法 ， 增 强 了 全 局 和 局 部 的 搜 
FCM 对 初始 聚 类 中 心 敏感 ， 需要 人 为 确定 聚 类 数 ， 容易 陷 入 局 。 索 能 力 和 效率 。 该 算法 在 收敛 时 间 略 有 增加 的 前 提 下 ， 可 以 有 
部 最 优 解 ， 同 时 受 噪声 点 及 孤立 点 对 聚 类 结果 的 影响 ， 使 混合 。 效 、 快 速 地 解决 FCM 算法 过 于 依赖 初始 聚 类 中 心 、 易 陷入 


和 


可 
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部 最 优 的 问题 ， 提 高 了 混合 矩阵 估计 结果 的 稳定 性 。 但 是 不 能 
克服 孤立 点 对 聚 类 中 心 的 影响 ， 为 此 本 文 又 利用 霍 夫 变换 来 修 


正 GASA-FCM 得 到 的 聚 类 结果 ， 有 效 地 提高 了 混合 矩阵 估计 
的 精度 。 
1 ”基于 信号 稀 疏 性 的 欠 定 混 合 矩 阵 估 计 
一 般 的 线性 混合 欠 定 盲 源 分 离 的 数学 模型 可 写 为 
x(1)=As(?) )=- Dasl 3 (1) 
其 中 : x(1)=[n (0), (0),.. ,0 表示 M 个 观测 信号 


(中 | 表示 N 个 源 信 号 ; 4e Rr 为 混 


s(7) =| 5,(7),s, (0),...,s 


稀 政 时 ， 信 号 在 绝 大 多 数 采样 点 处 幅 值 为 零 
同一 采样 时 刻 ， 只 有 一 个 源 信号 的 取 值 不 为 


合 矩 阵 。 当 源 信号 
或 者 接近 于 零 加 。 
零 即 


的 _ 太 (由 

uh ti 
但 在 实际 应 用 中 大 部 分 信号 在 时 域 中 并 不 是 充分 稀疏 的 ， 
因此 必须 寻求 其 他 变换 方法 增强 信号 的 稀 玻 性 四。 本 文采 用 短 


时 傅 里 时 变换 〈STFT) , 则 式 (1) 的 STEFT 表达 式 为 


AG) 


Ga; 


i 


xu (1) 
a 中 9 . 


X(t,m)=AS(t,w) ->as( 1,0@) (3) 


其 中 :XX(4,w) 和 ss,(1,w) 分 别 表示 x(1) 和 s(t) 的 STFT 变换 。 


党 号 经 过 STFT 变换 后 ， 式 〈2) 可 改写 为 


(aa EC X90) 00) 0 
CT (25 Gi Ca 
2 基于 GASAFCM 的 混合 聚 类 算法 
2.1 模糊 C- 均 值 聚 类 
模糊 C- 均 值 聚 类 属于 基于 目标 函数 的 模糊 聚 类 算法 的 范 
暑 中 。FCM 是 一 种 优化 算法 ， 它 通过 优化 目标 函数 求 出 每 个 数 
据点 隶属 度 ， 达 到 数据 点 分 类 的 结果 ，FCM 算法 简单 ， 且 精确 


度 较 高 ， 是 一 种 经 典 的 混合 矩 
假设 样本 集合 为 X = ( 划 ,x2,.….,xx) ,将 其 分 成 c 个 模糊 组 ， 

求 每 组 的 聚 类 中 心 c(i=1,2,.….,C)， 使 目标 函数 达到 最 小 。 
FCM 有 一 个 归 一 化 的 规定 , 即 整个 数据 集 的 隶属 度 总 和 为 


阵 估计 聚 类 算法 四。 


1。 
uj =1, 中 =1 (5) 
FCM 的 目标 函数 : 
TV.0..n0) = = -=> De (6) 
为 了 约束 目标 函数 ， 引 入 拉 格 朗 日 乘 子 和 ， 构 造 新 的 目标 
函数 
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J (Uc Ces hy hd )=T(U,c,..., Ce )+ D1 
n . [9 (7) 
=Y Yurd?+ Sm- 
i=l j j=: i=l 
对 上 式 求 偏 导 ， 得 到 目标 函数 的 最 小 值 ， 得 到 聚 类 中 心 : 
2 
Ci = 上 (8) 
2 
j=1 
隶属 度 : 
有 1 
i ad, 2/(m-1) (9) 
k=1 dy 
FCM 流程 如 图 1 所 示 
始 
vy 
参数 设置 ， 初 始 化 聚 类 中 心 
MU 
| 计算 初始 的 距离 矩阵 D 
e 
更 新 隶属 度 矩 阵 
机 
更 新 聚 类 中 心 
Vv 
重新 计算 距离 矩阵 ， 并 计算 目标 函数 的 什 


1 FCM 流程 


2.2 ”遗传 模拟 退火 算法 

遗传 模拟 退火 算法 是 一 种 基于 组 合 思想 而 产生 的 混合 优 
算法 钙 。 在 把 握 搜索 过 程 总 体能 力 比 较 强 的 遗传 算法 中 融入 
部 搜索 能 力 较 强 的 模拟 退火 算法 ， 增 强 了 全 局 和 局 部 的 搜索 
力 和 效率 。 这 两 种 算法 结合 后 ， 在 搜索 能 力 上 相互 补充 ， 是 一 


dl cc 


Tm 
CC 


种 在 性 能 上 更 优良 的 搜索 算法 ， 提 高 了 整个 算法 更 方面 的 搜 
能 力 ， 降 低 了 对 参数 的 要 求 n0， 有 具有 较 明 显 的 优越 性 。 
遗传 模拟 退火 算法 的 总 体 流 程 与 基本 的 遗传 算法 相 类 似 


系 


过 


先 产生 随机 的 初始 解 ， 利 用 遗传 算法 进行 全 局 搜索 ， 然 后 通过 
交叉 变异 等 操作 产生 新 个 体 ， 再 独立 地 对 产生 的 各 个 个 体 进 


模拟 退火 操作 , 得 到 下 一 代 种 群 个 体 。 重复 此 过 程 , 直至 满足 。 


基本 流程 如 下 
a) 初始 化 参数 〈 最 大 人 欠 代 次 数 、 遗 传 算法 交叉 率 、 模 拟 退 

火 算法 初始 温度 等 ); 

b) 随机 产生 初始 种 群 ; 


这 


何 . 
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c) J 


j 适 应 度 函 数 计算 初始 种 群 中 所 有 个 体 的 适应 度 ; 


d) 选择 种 群 中 的 个 体 进 
e) 将 模拟 退火 操作 运用 了 


行 交 叉 变 异 操作 ， 产 生 新 种 群 ; 


新 种 群 的 各 个 个 体 ,产生 新 个 体 ; 


f) 计算 个 体 的 适应 度 ， 根 据 适应 度 值 进行 评价 ; 
g) 判定 终止 条 件 ， 若 满足 则 输出 最 优 解 ， 算 法 结束 ; 若 不 


满足 则 转 到 步骤 d)。 


2.3 基于 SAGA-FCM 的 混合 聚 类 算法 


为 了 有 效 的 改善 FCM 


局 部 收敛 性 , 采 


] SAGA-FCM 的 混 


合 聚 类 算法 ， 通 过 模拟 退火 算法 较 强 的 局 部 搜索 能 力 和 遗传 算 


图 3 所 示 。 


开始 


法 较 强 的 全 局 搜索 能 力 ， 提 高 其 聚 类 结果 的 鲁 棒 性 。 其 流程 如 


初始 化 参数 ， 


生成 初始 种 群 


E 
vy 


计算 各 个 体 隶 


属 度 、 适 应 度 值 


P| 


v 


Gen=0 


一 一 


2 
遗传 算法 进行 选择 、 交 又 、 变 异 操作 


v 
新 和 


中 群 


聚 类 中 心 ， 


M4 
计算 各 个 体 的 隶属 度 、 适 应 度 值 和 c 个 


模拟 退 


火 算法 蔡 换 上 昌 个 体 


< 


T< 


N 


To 


N 


图 2 GASA-FCM 流程 


3 ” 霍 夫 变换 修正 聚 类 中 心 


霍 夫 变换 ， 图 形 空间 中 直线 


霍 夫 变换 是 将 图 形 空间 中 的 数据 变换 到 参数 空间 器。 经 过 
上 的 点 在 参数 空间 中 对 应 的 


线 相 


交 于 一 点 ， 交 点 处 会 形成 极 值 点 。 利 / 
的 主要 思路 是 : 将 时 频 域 中 的 观测 信号 进行 霍 夫 变换 ， 
上 的 观测 点 变换 得 到 的 所 有 


聚 类 直线 


j 霍 夫 变换 估计 混合 矩阵 
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上 六 es 
,Pp 
gE 
' ~ SS 
晤 
t 四 0 ] 2 
型 和 
[SEE [W343 
图 3 霍 夫 变换 原理 


主要 集中 在 n 条 直线 上 ， 直 线 向 量 : 


sng sn sin@, sin Oy ， 
cosO sin® sin@, sin Oy ， 
cos@ sin sin Oy _» 
coSsO ， 


式 (12) 中 性 =1 ，(8.2，… 


换 到 参数 空间 为 


=1'x 


过 


30) es[0， x)” 


于 语音 信号 存在 稀 琉 性 , 对 于 欠 定 情况 下 的 观测 数据 点 


sin Ov | 

sin Ow ， 

sin Oy ， GD 

sin Oy 

cosO | 

”。 将 直线 转 
(12) 


将 聚 类 得 到 的 结果 进行 霍 夫 变 换 修 正 聚 类 中 心 ， 弥 补 了 孤 


I 


立 点 对 估计 精度 的 影响 。 


混合 矩阵 估计 ， 首 先 对 时 频 域 的 观测 


本 文 提出 的 基于 GASA-FCM 混合 聚 类 和 霍 夫 变换 的 欠 定 
言 号 利用 GASA-FCM 进 
于 利用 霍 夫 变换 对 每 类 数 
十 计 值 。 本 文 提出 算 


行 聚 类 ， 得 到 混合 矩阵 估计 的 初 值 ; 
据 的 聚 类 中 心 进行 修正 ， 得 到 混合 矩阵 的 
法 的 具体 步骤 如 下 : 
a) 观测 信号 经 过 短 时 傅 里 叶 变 换 后 ,观察 时 域 和 频 域 的 散 
点 图 。 


b) 若 频 域 散 点 图 中 直线 方向 有 孤立 的 点 ， 采 


文献 [13] 提 


出 的 单 源 时 频 点 检测 方法 增强 信号 的 稀 玻 性 


c) 设 定 阔 值 去 除 低能 量 点 ,将 剩余 的 观测 信号 按 列 归 一 化 


同 
线 会 相交 形成 极 值 点 ， 


投影 到 M 维 上 半空 间 ， 观 察 观 测 
d) 采 


言 写 的 聚 类 性 


通过 对 极 值 点 进行 累加 得 到 局 部 极 大 值 点 的 位 
量 ， 进 而 得 到 混合 矩阵 的 估计 。 
假设 图 形 空间 里 过 点 (x,y) 的 直线 为 y=ax+b， 对 直线 进 


EE， 计 算 直线 向 


行 转换 到 变换 空间 p -9 里 的 参数 方程 如 下 : 


D=xcos(9)+ysin(O) 
其 中 : P 表 示 直 角 坐 标 系 中 原点 到 
直线 法 线 的 夹 角 。 霍 夫 变 换 的 原理 如 图 4 所 示 。 霍 夫 变换 将 图 
形 空间 里 同一 直线 上 的 点 A、B、C、D 变换 为 参数 空间 里 不 同 
线 A “、B'、C”、D' 且 相交 于 一 点 P”， 通 过 检测 P" 坐 


的 


参数 从 


标 位 


(10) 
直线 的 距离 ;0 表示 x 轴 与 


而 估计 直线 所 在 方向 [9。 


GASA-FCM 聚 类 算法 对 归 一 化 处 理 后 的 观测 信号 
进行 聚 类 ， 得 到 聚 类 中 心 和 初步 的 混合 矩阵 估计 值 


e) 对 经 过 GASA-FCM 聚 类 所 得 到 每 一 类 数据 进行 霍 夫 变 


换 ， 检 测 局 部 极 大 值 点 的 坐标 参数 ， 计 算 直线 向 量 ( 即 混合 矩阵 


列 向 量 )， 进 而 得 到 修正 后 的 混合 矩阵 估计 值 。 
4 ”实验 仿真 与 分 析 
4.1 混合 矩阵 评价 准则 


为 了 评价 欠 定 混合 矩阵 估计 的 精度 ， 本 文采 ) 


‘= 


归 一 化 均 方 误差 (normalized mean square error， 


文献 [14] 的 
NMSE) 和 文献 


[15] 的 偏离 角度 两 种 指标 来 衡量 所 提 算 法 的 估计 性 能 。 
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归 一 化 均 方 误差 表达 式 为 
2 (0a) 
NMSE =101ogi,| 之 己 (13) 


其 中 ;混合 矩阵 的 行 数 和 列 数 分 别 为 M、N; 混合 矩阵 4 和 估 
计 和 矩阵 A 第 ; 行 j 列 元 素 分 别 为 w 、a, 。NMSE 越 小, 则 估计 
的 混合 矩阵 越 精确 。 

偏离 角度 表达 式 为 


~、 180 < 0Gi > 
ang(a,, a, )= 0 E | | 2 | 


式 (14) 表示 wu 与 a 的 夹 角 。 其 中 : a 为 初始 矩阵 A 的 第 
i 列 ; 5 为 估计 矩阵 和 的 第 ; 列 (i=1,2,…,N )。 偏 离 角度 越 小 ， 
混合 矩阵 入 的 性 能 越 好 。 
4.2 实验 仿真 与 分 析 
为 了 验证 上 述 算法 的 有 效 性 ， 本 实验 在 MATLAB R2010b 
环境 下 ， 进 行 一 系列 仿真 对 比 实验 ， 源 信号 采用 NOIZEUS 语 
音 库 的 语音 信和 号， 采样 频率 为 8 kHz。 
仿真 1 将 三 路 源 信 号 混合 成 两 路 观测 信号 , 即 N=3,M=2。 
随机 生成 2X3 的 混合 矩阵 为 
nn 0.6207 
0.9317 -0.7841 


一 0.0848 
0.7288 


三 路 源 信号 经 过 混合 矩阵 A， 得 到 的 两 路 观测 信号 时 域 波 
形 如 图 4 所 示 。 对 观测 信号 进行 短 时 傅 里 叶 变换 , 选用 Hanning 
窗 ， 窗 口 长 度 为 512, 闭合 长 度 为 256。 变换 后 的 观测 信号 时 频 
或 散 点 图 如 图 5 所 示 。 


图 5 时 频 域 散 点 图 


| 
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图 5 可 知 ， 观 测 信 号 在 频 域 的 方向 性 不 清晰 ， 本 文采 用 
提出 的 时 频 单 源 点 检测 方法 增强 观测 信号 的 稀疏 性 ， 使 其 呈 3 
出 良好 的 方向 性 , 设 定 阔 值 == 0.15 去 除 低能 量 点 09。 对 剩余 
观测 信号 进行 聚 类 ， 得 到 混合 矩阵 的 估计 。 时 频 单 源 点 检测 后 
的 散 点 图 和 归 一 化 散 点 图 如 图 6 所 示 。 
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图 6 时 频 单 源 点 检测 后 的 散 点 图 和 归 一 化 散 点 图 


实验 参数 设置 为 暴 指 数 p=3; 最 大 达 代 次 数 T ”=100 ; 目 

标 函数 的 终止 容 限 o =1e-6; 冷 却 系数 q=0.6; 初 始 温度 到 =100; 

终止 温度 T=2; 最 大 遗传 代数 G， =10 。 为 了 说 明 本 文 算法 

的 有 效 性 ， 对 比 FCM、GAFCM、GASAFCM 和 本 文 算法 分 别 

所 得 到 的 混合 和 矩阵 。 
FCM 算法 估计 的 混合 矩阵 : 

-0 i 0.6159 | 


0.9266 -0.7872 0.7391 
GAFCM 算法 估计 的 混合 矩阵 : 


本 _「-0.3725 0.6220 -0.6752 
cc | 0.9280 -0.7830 0.7376 


GASAFCM 算法 估计 的 混合 矩阵 : 


品 —0.3654 0.6189 
AGasarcm = 0.9308 


—0.6829 
-0.7855 0.7305 


本 文 算 法 估计 的 混合 矩阵 : 
S bp 0.6213 


—0.6836 
0.9309 -0.7836 0.7298 


采用 归 一 化 均 方 误差 对 比 混合 矩阵 估计 的 精度 ， 分 析 算 法 
的 运行 时 间 ， 如 表 1 所 示 。 
表 1 各 算法 的 归 


化 均 方 误差 和 运行 时 间 的 对 比 


对 比 算法 ”NMSE/dB ”运行 时 间 /s 
FCM -38.1792 1.2237 
GAFCM -40.4446 1.4626 
GASAFCM -52.4875 2.2719 
本 文 算法 -56.1383 2.9528 


通过 计算 估计 混合 算 阵 与 初始 矩阵 每 一 列 的 偏离 角度 ， 进 
一 步 分 析 混 合 矩 阵 的 估计 精度 ， 如 表 2 所 示 。 

仿真 2 将 四 路 源 信号 混合 成 三 路 观测 信号 , 即 N=4,M=3。 
随机 生成 3X4 的 混合 矩阵 为 
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0.9543 -0.9748 
A=|-0.2220 -0.2149 0.6467 -0.7505 
0.2003 -0.0606 0.7342 -0.6318 


-0.2070 0.1936 


表 2 估计 的 入 与 初始 A 的 偏离 角度 对 比 /。 


a 与 的 偏离 角度 


对 比 算法 。 ang(a,d) 


ang(ayt) ang(a,t) 


FCM 0.7894 
GAFCM 0.5681 
GASAFCM 0.1308 
本 文 算法 0.1181 


0.3260 0.8760 
0.0975 0.7461 
0.1306 0.1460 
0.0447 0.0893 


| 


吧 


8 ”时 频 域 散 点 


9 时 频 单 源 点 测 后 的 散 点 图 和 归 一 化 散 点 


FCM 算法 估计 的 混合 矩阵 : 


GAFCM 算法 估计 的 混合 入 


| 0.9330 -0.9755 
Au =|-0.2171 -0.2148 
| 0.2447 -0.0473 


E 阵 : 


图 7 三 路 观测 信号 


尘 


对 


-0.1056 0.1477 
0.6782 -0.7505 
0.7269 -0.6078 


李 


ChinaXiv 合 作 期 刊 
虎 ， 等 : 基于 GASA-FCM 混合 聚 类 与 霍 夫 变 换 的 欠 定 混合 矩阵 估计 


0.9520 -0.9766 -0.1967 0.1821 
Are =| -0.2163 -0.2080 0.6583 -0.7372 
0.2166 -0.0545 0.7267 -0.6507 


GASAFCM 算法 估计 的 混合 矩阵: 


0.9570 -0.9782 -0.2094 0.1941 
A_wrey =|-0.2160 -0.1984 0.6444 -0.7530 
0.1938 -0.0616 0.7355 -0.6287 


本 文 算法 估计 的 的 混合 矩阵 : 


0.9533 -0.9757 -0.2118 0.1967 
A=|-0.2249 -0.2108 0.6458 -0.7500 
0.2016 -0.0640 0.7336 -0.6315 


表 3 各 算法 的 归 一 化 均 方 误差 和 运行 时 间 的 对 比 
对 比 算法 ”NMSE/dB ”运行 时 间 /s 


FCM -23.8101 1.4990 
GAFCM -34.7021 1.7417 
GASAFCM -40.0089 2.6494 
本 文 算法 -48.0329 3.4143 


表 4 估计 的 入 与 初始 A 的 偏离 角度 对 比 /。 


a 与 亏 的 偏离 角度 


对 比 算法 ang (a,é) ang (a,,t,) ang (as,6,) ang (as,,) 
FCM 2.7756 0.7631 6.1014 2.6912 
GAFCM 0.9981 0.5376 0.9872 1.4791 
GASAFCM 0.5299 0.9669 0.2045 0.2300 
本 文 算法 0.1909 0.2408 0.2819 0.1807 


表 1 和 3 混合 矩阵 的 NMSE 和 运行 时 间 可 知 , 本 文 提出 


的 算法 虽然 NMSE 比 其 他 几 种 算法 小 , 有 旦 混合 矩阵 估计 的 精度 


越 准 确 ， 但 算法 的 运行 时 间 明 显 比 其 他 几 种 算法 长 。 而 由 表 2 


和 4 


可 知 ， 本 文 算法 得 到 的 混合 矩阵 每 列 的 偏离 角度 均 最 小 ， 
估计 精度 更 准确 ， 与 表 1 和 3 分 析 所 得 结论 一 致 ， 本 文 算法 能 
够 有 效 、 准 确 地 估计 混合 矩阵 。 


为 了 验证 本 文 算法 的 实用 性 ， 通 过 大 量 实验 数据 及 仿真 1 


和 2， 


分 析 NMSE 对 比如 图 10 所 示 。 


归 一 化 均 方 误差 
ww 二 


-60 —A—FCM 
—E—GAFCM 
—e— GASAFCM 
一 一 本 文 佣 法 

3 4 5 6 


10 NMSE 对 比 


图 10 可 知 ， 随 着 源 信号 个 数 的 增加 ， 估 计 的 混合 矩阵 


NMSE 不 断 增 大 。 
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5 ”结束 语 


sn ner ii ee 针 
对 传统 的 FCM 算法 过 于 依赖 初始 聚 类 中 心 ， 搜 索 过 程 易 陷入 


局 部 最 优 ， 本 文 将 遗传 算法 与 模拟 退火 算法 相 结合 ， 互 相 取 长 
补 短 ， 提 出 一 种 基于 GASA-FCM 混合 聚 类 与 霍 夫 变换 的 欠 定 


混合 矩阵 估计 算法 ， 增 强 了 局 部 和 全 局 搜索 能 力 ， 有 效 改善 
传统 FCM 算法 容易 陷入 局 部 最 优 的 缺陷 ， 同 时 ， 为 了 提高 混 
合 矩 阵 估计 的 精度 ， 对 聚 类 结果 引入 霍 夫 变换 修正 聚 类 中 心 。 
仿真 实验 结果 表明 ,该 算法 提高 了 混合 矩阵 估计 精度 和 算法 
棒 性 ， 对 混合 矩阵 的 估计 是 可 行 和 有 效 的 。 
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